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By means of the method of the effective potential, the phase transitions from the Mott insulating 
state to either the atomic or the dimer superfiuid state in the three-body constrained attractive Bose 
lattice gas are analyzed. Due to the appearance of the Feshbach resonance coupling between the two 
kinds of order parameters in the derived effective potential function, it is found that the continuous 
Mott insulator-to-superfluid transitions can be preempted by first-order ones. Since the employed 
approach can provide accurate predictions of phase boundaries in the strong coupling limit, where 
the dimer superfiuid phase can emerge, our work hence sheds light on the search of this novel phase 
in real ultracold Bose gases in optical lattices. 
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The impressive developments on the manipulation of 
ultracold gases in optical lattices provides one of the best 
environments in the search for exotic quantum phases [l| . 
Because unprecedented control over microscopic model 
parameters can be achieved, it is possible to explore pa- 
rameter regimes which are not available in other analo- 
gous condensed matter systems. The remarkable experi- 
mental demonstration of the superfiuid to Mott insulator 
transition in ultracold lattice bosons has paved the 
way for investigating other strongly correlated phases in 
similar setups. For instance, the search for novel and 
unconventional quantum phases in the mixtures of ultra- 
cold atoms obeying either the same or different statistics 
has attracted considerable attention 

It was recently suggested that intriguing quantum crit- 
ical behaviors can occur in attractive bosonic lattice gases 
with three-body on-site constraint [1, 3] . The system is 
described by the Bose-Hubbard model with a three-body 
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Here, ai{a\) is the bosonic annihilation (creation) opera- 
tor at site i, t IS the hopping matrix element, f/ < the 
on-site two-body attraction, and /i the chemical poten- 
tial. The convention signifies a sum over nearest- 
neighbor sites i and j. The on-site constraint can arise 
naturally due to large three-body loss processes d, Q, 
and it stabilizes the attractive bosonic system against 
collapse. Therefore, besides the conventional atomic su- 
perfiuid state (ASF) with non-vanishing order parame- 
ters (a) ^ and (a^) ^ appearing in the weakly- 
interacting limit, a dimer superfiuid phase (DSF) with 
vanishing atomic order parameter ((a) = 0) but nonzero 
pairing correlation {{a?) ^ 0) can be realized for suffi- 
ciently strong attraction Q. It was shown in Refs. 0, Q 
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that this model provides a simple realization of the 
physics of Ising quantum transition together with the 
Coleman- Weinberg mechanism Q without resorting to 
the Feshbach-resonant mechanism While the na- 

ture around the ASF-DSF transition has been discussed, 
the detailed physics of the Mott insulator (MI) to super- 
fiuid (either ASF or DSF) transitions is not addressed in 
Refs. 0,11. 

In the present work, we focus our attention on the 
MI- ASF and the MI-DSF transitions in this three-body 
constrained attractive Bose lattice gas. Since there is no 
hopping term for dimers in the model of Eq. ([T}, usual 
strong-coupling theory based on simple mean-field decou- 
pling is not appropriate for our purpose, because it fails 
to describe the DSF phase. Instead, the method of the ef- 
fective potential developed in Refs. [ll|, [l^l is employed, 
which has been applied successfully to the usual Bose- 
Hubbard model with repulsive interaction. Accurate ana- 
lytical results for the phase boundaries of Mott insulator- 
to-superfluid transitions have been obtained when the 
contributions from higher orders in the hopping term are 
included systematically [ll,[l3l- Here, we generalize their 
approach to the case with two order parameters. For the 
attractive bosonic lattice gases under consideration, we 
show that the effective potential function has the same 
form as the mean-field Ginzburg-Landau theory of reso- 
nant Bose gases [8l-[l0|. Following similar discussions in 
Ref. [l^l, it is found that the MI-DSF transitions are al- 
ways second-order, while the MI-ASF transitions can be 
either second-order or first-order. The tricritical points 
on the MI-ASF transition lines and the critical end points 
of the MI-DSF phase boundaries are determined within 
the present approach. From our results, it is shown that 
the DSF phase exists only in a narrow region of chemi- 
cal potential ^/\U\ for small hopping parameters t/\U\. 
Since the approach used in the present work can be con- 
sidered as a kind of strong-coupling expansion, it is ex- 
pected to provide accurate phase boundaries of the MI 
phases for the transitions occurring at small values of 
t/\U\. Our results hence provides a useful guide to the 
experimental search of the DSF phase and the associ- 
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ated quantum phase transitions in ultracold Bose gases 
in optical lattices. 

To derive the effective potential with two order pa- 
rameters for the model in Eq. (IlJ, we begin by adding 
two symmetry breaking source terms to the Hamiltonian, 
which are spatially and temporally global: H ^ Hq + V 
with Hq being the single-site zero-hopping contribution 
in Eq. © and T/ = a^j + J:^{x*a^ + r^^af + 

H.c). It is easy to see that the exact ground state for 
the on-site part Hq is the n = MI state (i.e., the empty 
state) for the chemical potential /i < —\U\/2, while it 
becomes the n = 2 MI state (i.e., the completely filled 
state) for /i > — 1[/|/2. By treating F as a perturbation 
and following the specific adaption of high-order many- 
body perturbation theory proposed in Refs. [HI [13: the 
free energy (or the ground state energy at the T = limit 
under consideration) in the MI phases for the modified 
Hamiltonian H can be calculated as a double power se- 
ries in both the hopping parameter t and the source fields 
X, X* 3.nd rj, 77*. Up to the forth power in source fields, 
the general expression of free energy per site / takes the 
following form 



while, for the n = 2 MI state with /i > —1/2, 
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where /o corresponds to the ground state energy density 
in the absence of the source fields. We point out that 
our expression of / does respect the symmetry for the 
U(l) phase transformation X X e*^ ^nd 7y — > 77 e^**, 
which can be understood from the form of the modified 
Hamiltonian H . 

The coefficients in the expansion of free energy density 
/ in Eq. ^ can be determined perturbatively in hopping 
parameter t. In the following, energy unit is set to be \U\ 
and we define the dimensionless parameters t = t/\U\ 
and pL — ^Ji./\U\ for convenience. Up to leading order in t, 
we find that, for the n~Q MI state with jl < —1/2, 
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Here z is the coordination number of the underlying lat- 
tices. Form the free energy density /, the order parame- 
ters of the atomic condensate (p^. = (a) and the molecular 
(or dimer) condensate 4>-m = (a^) can be obtained by the 
first derivative of / with respect to corresponding exter- 
nal sources. That is, 0a = df /dx* and 0m = df/d?f, 
respectively. We note that, due to the mixing terms in 
/ with coefficients A and U12, there exists a nontrivial 
relation between the two order parameters 0a and 0ni- 

The effective potential in terms of these order parame- 
ters is then derived by performing the Legendre transfor- 
mation on the free energy density /, r(0a, 4>l.^ 0m, 0m) = 
f ~ X*4>a. - X4>1 - '7*0m - '^0mi which Can be used to 
determine the phase boundaries of the Mott insulator- 
to-superfluid transition and their nature as shown below. 
From Eq. ^ , the Ginzburg-Landau expansion of the ef- 
fective potential as a power series of the order parameter 
variables can be obtained: 



''a, H'n, Vm, H^n-L 



h + mi|0aP + ||0a|' + m2|0mP + ^|0m|" 



-(0:;0^+C.C.)+5l2|0an0„.P 



(5) 



with the coefficients given by mi = — l/ri, gi = ui/rf + 
3X^/2rfr2, m2 = -l/r2, 92 = W2/?'2j " = X/rfr2, and 
512 = SX"^ /rir2—3ui2/rfrl. An interesting feature of the 
two-component Landau theory in Eq. ^ is the appear- 
ance of the Feshbach resonance coupling with coefficient 
a. The existence of such a term has important conse- 
quences on the Ising quantum phase transition in the 
boson Feshbach resonance problems [1, |^ . Starting from 
a specific mean-field ansatz for the DSF state, the anal- 
ogy between the present problem and the usual boson 
Feshbach resonance model was also noticed and explored 
in some detail in Ref. Q. Here we approach this issue 
from the MI states, and focus our attention on its impli- 
cation on the MI-ASF and the MI-DSF transitions. 

As usual discussions for Ginzburg-Landau theory, the 
continuous MI-DSF and MI-ASF transitions occur when 
the coefficients of the quadratic terms for the correspond- 
ing order parameters vanish. That is, 1712 = and 



3 



mi = give continuous MI-DSF and MI-ASF transi- 
tions, respectively [l^. From the relation between 7712 
and r2 derived above and the perturbative results for 
r2 in Eqs. ([3]) and the critical value of hopping pa- 
rameter for the continuous n = MI-D SF transition is 
given by 2;tc,oMi-DSF — \^ zp.{2fl + l)/2, while that for 
the continuous n = 2 MI-DSF transition is z ic,2Mi-DSF — 
y/z{jl + l)(2/2 -I- l)/2. Similar reasoning applies also for 
the continuous MI- ASF transition. Our calculations give 
the critical values of t for the continuous n = MI- ASF 
transition: zt^fiui-ASF — smd that for the continu- 
ous n ~ 2 MI- ASF transition: 2;tc.2Mi-ASF — (m + l)/2. 
The phase boundaries of continuous MI-DSF and MI- 
ASF transitions for bosons hopping on two-dimensional 
(2D) square lattices with coordination number z = 4 are 
depicted in Fig. [Ija) as the red and the blue solid lines, 
respectively. 

As mentioned before, the presence of the Feshbach res- 
onance term plays an important role in characterizing 
the nature of the ASF state, and it may even modify the 
feature of the phase transitions. In general, the phase 
boundaries of the MI phases are determined by minimiz- 
ing the effective potential F in Eq. ([5]) , which leads to the 
following extremum equations: 
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From Eq. (j6bp . it is clear that, once 0a becomes nonzero, 
the molecular condensate (pra will acquire a non- vanishing 
value also. Thus the fact that the atomic condensate 
and the molecular condensate are both finite in the ASF 
phase can be explained straightforwardly as long as the 
Feshbach resonance term exists in F. Nevertheless, since 
(/>a = in both the DSF and the MI phases, one can 
realize again from Eq. (j6bp that the condition for the 
continuous MI-DSF transition is not affected by the pres- 
ence of this Feshbach-resonancc term. Similarly, far away 
from the ASF-DSF phase boundary, 0m ~ in the ASF 
phase and = in the MI phase, thus the effect of the 
Feshbach resonance coupling on the aforementioned con- 
tinuous MI-ASF transition should be negligible, which 
can be understood from Eq. ([6a|) . 

However, around the parameter region such that both 
mi and m2 are zero, the continuous Ml-to-superfluid 
transitions mentioned above will be preempted by first- 
order MI-ASF transitions caused by the Feshbach res- 
onance term. Therefore, crossing this phase boundary, 
the atomic and the molecular condensates as well as the 
particle density will no longer change continuously. The 
occurrence of this first-order MI- ASF transition seems to 
be overlooked in Refs. 0, Q. As explained below, this 
conclusion can be reached by carefully analyzing the ex- 
tremum equations in Eq. (pa|) and ((6b| . which has been 
discussed in the context of the boson Feshbach resonance 
problem [Toj . 

Without loss of generality, we may take both 0a and 
0in to be real. We remind that, on the first-order MI- 
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FIG. 1: (Color online) Phase diagram of the three-body 
constrained attractive Bose-Hubbard model on (a) two- 
dimensional square lattices (b) three dimensional cubic lat- 
tices. The red solid lines indicate the continuous MI-DSF 
transitions, and the blue solid (dashed) lines show the continu- 
ous (first-order) MI- ASF transitions. The black dots (crosses) 
denote the tricritical (critical end) points. The schematic 
phase boundary between the ASF and the DSF phases is 
added as the green dashed line for clarity. The insets show 
the details of the first-order n = 2 MI-ASF transition lines. 



ASF phase boundary, there should exist a discrete jump 
for 0a to a finite value. Therefore, to determine this 
first-order phase boundary, we need to seek for non- 
trivial solutions of 0a. From Eq. (|6a[) . this is given by 
|0ap = (a0m — ,9i2|0mp — m-i)/^!- Substituting this Con- 
dition into Eqs. ()6bp . one leads to a saddle-point equa- 
tion, which is expressed solely in terms of 0ni, 
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Besides, the corresponding energy density difference be- 
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tween the ASF and the MI phases is given by 
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The first-order MI- ASF phase boundaries are determined 
by the simultaneous solutions of Eq. ([7]) and AF = in 
Eq. (|8]). Moreover, the tricritical points separating the 
continuous and the first-order MI-ASF transitions can 
be found by satisfying the conditions for both kinds of 
transitions. That is, by requiring nontrivial solution of 
Eq. ([7]) with mi = 0, one leads to the condition for the 
tricritical points: m2 — a^/2gi = 0. From the perturba- 
tive results for the coefficients of the free energy density 
/ in Eqs. ([3]) and ^ as well as the relations between 
the coefficients in / and those in the effective potential 
F, these first-order transition lines and the correspond- 
ing tricritical points can be calculated numerically. Our 
results of the first-order MI-ASF phase boundaries for 
the models defined on 2D square lattices with z = 4 
are depicted as the blue dashed lines in Fig. [Ija). The 
value of the tricritical point on the n = MI-ASF tran- 
sition line is found to be {It, fir) — (0.226,-0.903), 
while (iT,MT) ^ (0.071,-0.430) for that on the n = 2 
MI-ASF transition line. On the other hand, the value 
of the critical end point of the n = MI-DSF transi- 
tion line is given by (^EiMe) — (0.138,-0.624), while 
{tE,jiE) ^ (0.061, -0.472) for that of the n = 2 MFDSF 
transition line. These tricritical (critical end) points arc 
denoted by black dots (crosses) in Fig.[TJa). 



Within the present approach, the dependence of 
the lattice structure and the dimensionality enter only 
through the coordination number z. As an illustration 
for its influence, the results for the models on three- 
dimensional cubic lattice with z = 6 are presented in 
Fig.[TJb). We find that the ranges in p. (and in t also) for 
both the first-order transition lines and the DSF phase 
arc narrower than those in the 2D case [l^ . This implies 
that the role played by the Feshbach resonance term be- 
comes less important as z increases. 

To summarize, the phase boundaries of the MI-ASF 
and the MI-DSF transitions in the three-body con- 
strained attractive Bose lattice gas are determined by 
generali zing the effective potential approach developed 
in Refs. [ll[[l3]- While the coefRcients of the free energy 
density / are calculated only up to leading order in t/\U\, 
our results should be quantitatively accurate at least in 
the strong-coupling region of the phase diagram. We find 
that, due to the presence of the Feshbach resonance term 
in the Ginzburg-Landau expansion, the continuous MI- 
to-superfluid transitions can be driven to be first-order 
ones. From our results, it is found that the DSF phase 
exists only in a narrow region of chemical potential fJ./\U\ 
for small hopping parameters t/\U\. Therefore, carefully 
tuning system parameters into the suggested parameter 
regime are necessary to uncover experimentally this novel 
phase in real ultracold Bose gas in optical lattices. 

We are grateful to K.-K. Ng for for many enlightening 
discussions. Y.-W. Lee and M. -F.Yang thank the support 
from the National Science Council of Taiwan under grant 
NSC 96-2112-M-029-006-MY3 and NSC 96-2112-M-029- 
004- MY3, respectively. 



[1] M. Lewenstein, A. Sanpera, V. Ahufinger, B. Damski, A. 
Sen, and U. Sen, Adv. Phys. 56, 243 (2007); I. Bloch, 
J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 
(2008). 

[2] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, and 

I. Bloch, Nature (London) 415, 39 (2002); T. Stoferle, H. 

Moritz, C. Schori, M. Kohl, and T. Esslinger, Phys. Rev. 

Lett. 92, 130403 (2004); I. B. Spielman, W. D. Phillips, 

and J. V. Porto, Phys. Rev. Lett. 98, 080404 (2007); ibid. 

100, 120402 (2008). 
[3] S. Diehl, M. Baranov, A. J. Daley, and P. ZoUer, Phys. 

Rev. Lett. 104, 165301 (2010). 
[4] S. Diehl, M. B aranov, A. J. Daley, and P. Zoller, 

larXiv:0912.3 T92. and 0912.3196. 
[5] A. J. Daley, J. Taylor, S. Diehl, M. Baranov, and P. 

Zoller, Phys. Rev. Lett. 102, 040402 (2009). 
[6] M. Roncaglia, M. Rizzi, and J. I. Cirac, Phys. Rev. Lett. 

104, 096803 (2010). 
[7] S. Coleman and E. Weinberg, Phys. Rev. D 7 1888 

(1973); B. I. Halperin, T. C. Lubensky, and S.-K. Ma, 

Phys. Rev. Lett. 32, 292 (1974). 
[8] L. Radzihovsky, J. I. Park, and P. B. Weichman, Phys. 

Rev. Lett. 92, 160402 (2004). 
[9] M. Romans, H. Stoof, and S. Sachdev, Phys. Rev. Lett. 



93, 020405 (2004). 
[10] K. Sengupta and N. Dupuis, Europhys. Lett. 70, 586 

(2005); L. Radzihovsky, P. B. Weichman, and J. I. Park, 

Ann. Phys. 323, 2376 (2008). 
[11] F. E. A. dos Santos and A. Pelster, Phys. Rev. A 79, 

013614 (2009). 

[12] B. Bradlyn, F. E. A. dos Santos, and A. Pelster Phys. 
Rev. A 79, 013615 (2009). 

[13] N. Teichmann, D. Hinrichs, M. Holthaus, and A. Eckardt, 
Phys. Rev. B 79, 100503(R) (2009); ibid. 224515 (2009). 

[14] A. Eckardt, Phys. Rev. B 79, 195131 (2009). 

[15] The part of continuous Ml-to-superfluid transitions in 
the present model should belong to the z = 2 dilute bose 
gas transition universality class. See , for example, S. 
Sachdev, Quantum Phase Transitions (Cambridge Uni- 
versity Press, Cambridge, UK, 1999). 

[16] Now the value of the tricritical point on the n = MI- 
ASF transition line becomes {iT,fiT) — (0.119,-0.713), 
while {iT,flT) ^ (0.045,-0.456) for that on the n = 2 
MI-ASF transition line. On the other hand, the value of 
the critical end point of the n — MI-DSF transition line 
is given by (tE,/iE) — (0.083,-0.572), while (iE,A*E) — 
(0.039, -0.482) for that of the n = 2 MLDSF transition 
line. 



